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(Received ) 

The topological charge distribution P(Q) is calculated for lattice CP^ -1 models. In 
order to suppress lattice cut-off effects, we employ a fixed point (FP) action. Through 
transformation of P(Q), we calculate the free energy F(9) as a function of the 9 parameter. 
For N=4, scaling behavior is observed for P(Q) and F(0), as well as the correlation lengths 
£(Q). For N=2, however, scaling behavior is not observed, as expected. For comparison, 
we also make a calculation for the CP 3 model with a standard action. We furthermore pay 
special attention to the behavior of P(Q) in order to investigate the dynamics of instantons. 
For this purpose, we carefully consider the behavior of Jeff, which is an effective power of 
P(Q) (~ exp(— CQ 7ci? )), and reflects the local behavior of P{Q) as a function of Q. We 
study 7egt for two cases, the dilute gas approximation based on the Poisson distribution of 
instantons and the Debye-Hiickel approximation of instanton quarks. In both cases, we find 
behavior similar to that observed in numerical simulations. 



§1. Introduction 

It is interesting to study the phase structure of asymptotic free theories such 
as QCD and the CP^ -1 model. Non-perturbative studies of the phase structure 
of such theories are necessary in order to understand why effects of the topological 
term [9 term) are suppressed in Nature. The term affects the dynamics at low 
energy and is expected to lead to rich phase structures. ^ Actually, in the Z(N) gauge 
model, it has been shown by use of free energy arguments that oblique confinement 
phases could emerge and that an interesting phase structure may be realized. 2 ) In 
this paper we are concerned with the dynamics of the vacuum of CP^" 1 models 
with a topological term, which have several dynamical properties in common with 
QCD. We believe that study of the two-dimensional model will be useful in acquiring 
information about realistic physics. 

From the numerical point of view, the topological term introduces a complex 
Boltzmann weight in the Euclidean lattice path integral formalism. The complex 
nature of the weight prevents one from straightforwardly applying the standard 
algorithm used for Monte Carlo simulations. This problem can be circumvented 
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by Fourier-transforming the topological charge distribution P(Q). 3 '' 4 '' It is then 
necessary to calculate P(Q) as precisely as possible in order to reduce errors in the 
expectation values of physical operators as functions of 9. The precise determination 
of P(Q) will allow us to obtain a partition function and other quantities as functions 
of 9 with high precision. 

It is non-trivial how lattice cut-off effects would emerge through the introduction 
of a topological term. In the present paper, we employ fixed point (FP) actions to 
study this issue. 5 ** In the case of no topological term, the FP action is known to 
significantly suppress lattice cut-off effects for topological objects in CP^ -1 mod- 
els. 6 ) " 8 ) In Ref. 7 \ a FP action for the CP 1 model (CP^P) was determined, and the 
stability of instantons under minimization of the action was investigated in detail. 
It was also observed there that dislocations can be eliminated by adopting a FP 
charge as well as the FP action. However, the scaling behavior of the lattice topo- 
logical susceptibility xt was found to be strongly violated even after the dislocations 
are eliminated. For the CP 3 model with an FP action (CP 3 FP), contrastingly, im- 
pressive improvements have been found; 8 ) after the topological defects are removed, 
clear scaling behavior of xt is observed. 

In the present paper, we study the topological terms of CP 3 and CP 1 models 
with a FP action. For comparison, we also carry out a calculation for the CP 3 
model with the standard action (CP 3 ST). The main purposes of the paper are the 
followings. 

1. To study the scaling behavior of various quantities such as P(Q) , the free energy, 
the expectation value of the topological charge, the topological susceptibility, 
and correlation length function of Q. 

2. To analyze P(Q) in detail by considering an effective power j e jj of lnP(Q). 

The second of these purposes is associated with the phase structure of the model. 
In the very strong coupling region, there exists a first-order transition at 9 = it. 9 ) 
In this region P(Q) is Gaussian, and its volume dependence is given by P(Q) = 
C exp(— a/VQ 2 ), where C and a are /^-dependent constants, (3 being the inverse 
coupling. This 1/V-law of the exponent is associated with the existence of the first- 
order phase transition at 9 = 7T. 10 )' 11 ) When (5 becomes larger for a fixed volume, 
P(Q) has been found to deviate from the Gaussian form. As a consequence, the 
singularity at 9 = ir is no longer visible. It is interesting to consider the fate of the 
first-order phase transition. In order to see how P(Q) varies, we define an effective 
power P(Q) oc exp(— CQ" /e ff ), which is defined using three adjacent values of Q in the 
whole range of Q. We investigate the behavior of J e ff{Q) by systematically varying 
(3 and V. For a fixed value of f3, j e ff(Q) exhibits interesting behavior as a function 
of the topological charge density Q/V. For a small lattice size L, j e ff{Q) approaches 
some asymptotic value from below, while for large L it does so from above. For the 
range investigated, j e ff{Q) is always bounded from below by 1. As long as the finite 
size effects are not significant, jeffiQ) is between 1 and 2. Finite size effects are 
clearly seen in the behavior of j e ff(Q) when it exceeds 2. 

The second purpose stated above is also associated with the dynamics of instan- 
tons. We attempt to extract some information about these dynamics from ^ e g. For 
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this purpose, it is useful to employ two analytical models. One is that of a dilute 
gas of instantons obeying the Poisson distribution. For values of the parameter cor- 
responding to the very strong coupling region, the Poisson distribution leads to a 
Gaussian form of P(Q), but P(Q) deviates from the Gaussian form as the coupling 
constant becomes weaker. The other model is the Debye-Hiickel (D-H) approxima- 
tion of an instanton quark gas. 12 ) This is based upon an instanton quark picture, 13 ) 
in which instanton quarks interact weakly with each other. For these two models, 
P(Q) is calculated from the partition function Z(N+,N—), which is the probability 
to generate N + instantons and N_ anti-instantons. The quantities "f e ff(Q) for the 
two models exhibit behavior similar to that found in Monte Carlo simulations. Our 
conclusion is that the distribution P(Q) generated by Monte Carlo simulations are 
not inconsistent with the dilute gas approximation. 

In the following section, we briefly summarize the notation and the algorithm of 
the complex action calculation. In section 3 we present the results for CP 3 FP, and 
in section 4 we compare them with those of CP 1 FP and CP 3 ST. In section 5, we 
summarize and discuss analytical results for the Debye-Hiickel approximation of an 
instanton quark gas. We also compare values of j e jj from numerical simulations with 
those obtained from the Debye-Hiickel model and from the Poisson distribution. A 
summary is given in section 6. 

§2. Formulation 

2.1. Definition and algorithm 

The action with the 6 term is defined by 

S e = S-i0Q, (2-1) 

where S is a lattice action of the CP 7 ^ 1 model. Among various definitions of the 
topological charge, we here choose the geometrical definition. 14 ) The topological 
charge Q is counted by A n as 

in the updating process, where the plaquette contribution An is given by 

A\j = — ^{A M (n) + A v (n + n) - A^n + u)-A u (n + fi)}e^ u mod 2tt. (2-3) 

Here exp(iA^(n)) = z^{n)z(n + n)/\z^(n)z(n + i.e., A^(n) = aig(z^(n)z(n + n)). 
The link variables A^{n) = arg(z^ (n)z(n + /x)) satisfy 

Afi(n) G [— 7r, it]. (2-4) 

In order to avoid the complex Boltzmann weight, we adopt an algorithm by which 
the partition function is given by the Fourier transform of the topological charge 
distribution P(Q) 

Z(8) =Y j e i6Q P{Q). (2-5) 

Q 
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The distribution P(Q) is given in terms of the real Boltzmann weight as 

P(Q) = j [dzdz} Q e~ s / J [dzdz]e- s , (2-6) 

where [dzdz]® is the constrained measure on which the value of the topological charge 
given in Eq. (2-2) is restricted to Q. P(Q) is normalized as J2qP(Q) = 1- 
The expectation value of an operator O as a function of 9 is given by 

(O(0)) = J2 e Mi0Q]P(Q)(O) Q /Z(9), (2-7) 
Q 

where 



(O)q = J [dzdz] Q Oe- s / J [dzdz}e- s . (2-8) 



We calculate P(Q) by updating configurations through the combined use of the over- 
relaxation algorithm and the Metropolis algorithm for CP 3 , while only the Metropo- 
lis algorithm is applied for CP 1 . From the generated configurations, the topological 
charge is calculated according to Eq. (2-2). Since each function P(Q) under con- 
sideration rapidly falls off, it is convenient to restrict the range of Q for a single 
Markoff chain. We use the set method 15 )' 16 ) by which an entire range of Q is di- 
vided into sets Si,i = 1,2,3, Typically, each of the sets Si consists of 4 bins, 
Q = 3i — 3, 3i — 2, 3i — 1 and 3i, so that the adjacent set overlaps at the edge bin 
of the set, Q = 3k (k = 1, 2, 3, ■ ■ •). Depending on (3, the volume and also N, it is 
sometimes more convenient to use a wider range of bins for a set to save computation 
time, as well as to allow for better use of the trial function method described below. 

In order to generate configurations more efficiently, an effective action is used. 16 ) 
This action is modified by adding a trial function Pt(Q) according to 

S eff = S-\nP t (Q). (2-9) 

The form of Pt(Q) is chosen to be 

P t (Q) oc exp(-aQ^), (2-10) 

where a and 7 are adjusted so that P(Q) becomes almost flat, in order to reduce 
errors. The power 7 is often chosen to be 2.0 (the Gaussian case). 



2.2. Fixed point action 

The idea of using a FP action in asymptotically free theories in order to remove 
lattice artifacts was introduced by Hasenfratz and Niedermayer. 5 ) The action is 
defined at the fixed point of a renormalization group transformation at (3 = 00 and 
is perfect, i.e. without any discretization errors, on the classical level. Although on 
the quantum level the FP action is not perfect, it is in practice powerful with respect 
to removing lattice defects. For a finite coupling constant (3, the action S is given 
by 

S = NpA FP , (2-11) 



CP Models with a 9 Term and Fixed Point Action 5 



where ^4fp is the FP action. App is determined by a block transformation at f3 = oo 
to be 

A FP (0=min z [A FP (z)+T(C,z)], (2-12) 

with the transformation kernel T for block transformation from CP^ -1 spins z to 
block spins £ (see Ref . 5 ) for details) . It has been shown that lattice defects are 
invisible up to fairly small coupling constants that correspond to correlation lengths 
of only a few units of the lattice spacing. Because of these benefits we use the FP 
action to study the CP 3 model with a topological term. 

For numerical simulations, a parameterized form of the FP action is required. In 
previous studies, a large number of coupling constants have been used for parametriza- 
tion. In Ref. 8 ) , for example, CP 3 FP has been parameterized using 32 coupling con- 
stants in order to realize high precision. In order to reduce the computational effort 
required to carry out simulations with a FP action, we constructed for the present 
work a simpler parametrization using the same method and the same set of con- 
figurations as in Ref. 8 ). We were able to reduce the number of coupling constants 
from 32 to 9, while only increasing the average relative deviation between the mini- 
mized and the parameterized FP action from 0.4% to 0.6%. All the couplings of the 
parametrization are limited to a short range and lie within one plaquette. We list 
the coupling constants of the simpler parametrization in Table I, where we use the 
same numbering scheme as in Ref. 8 ) . 

For CP 1 FP, we employ the same set of coupling constants as in Ref. 5 ) (24 types 
of coupling constants). 



Table I. Coupling constants of the CP FP action used in this work. 



No. in Ref. 8) 


Coupling 


No. in Ref. 8) 


Coupling 


No. in Ref. 8) 


Coupling 


1 


0.61884 


2 


-0.04879 


4 


0.19058 


7 


0.00645 


11 


0.02984 


17 


-0.12433 


19 


0.04034 


20 


0.14477 


23 


0.01852 



2.3. Measurements 

The free energy density F{9) is obtained from the partition function of Eq. (2-5) 
through the relation 

F(0) = ~]nZ(6), (2-13) 

where V = L 2 , and L is a dimensionless lattice extension. The expectation value of 
the topological charge is defined as 

«>. = -H)^- (2-14) 

The correlation length £(Q) = l/m(Q) for a fixed topological charge sector is 
obtained from two point functions of P = z <S> z projected to zero momentum. It is 
extracted by analyzing their long distance fall-off of the form 

(P(0)P(r)) Q = A (exp (-m(Q)r) + exp (-m(Q) (T - r))) , (2-15) 

where T is the extent of the lattice in the time direction. 
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§3. Numerical results for CP 3 FP 

3.1. CP 3 FP without a topological term 

Before presenting the results of our simulations using the set method, we briefly 
discuss numerical results of a series of standard simulations using the FP action for 
the CP 3 model. There are two reasons why it is useful to perform these additional 
simulations. First, it has been observed 8 ) that use of a FP action leads to good 
scaling behavior of the topological charge even with a coarse lattice. Since in the 
present work we use a simpler parametrization of the FP action, we have to check 
that these good scaling properties are kept for this new action. The second reason is 
that from simulations without a topological term basic information about correlation 
lengths is obtained. This is needed to answer fundamental questions such as which 
couplings should be considered as strong and which as weak and which values of the 
lattice size L correspond to a large lattice in physical units and which to a small one. 
Values of the correlation length are also required in order to allow for investigation of 
scaling behavior. In the following sections it is understood, if not otherwise stated, 
that the correlation length is that determined from simulations without a topological 
term. 

Table II. Results of simulations with CP FP without a topological term. 



p 


L 


a 




XST 


Xfp 




xW 2 


0.5 


10 


0.4515(61) 


22.15(30) 


0.05654(11) 


0.05394(11) 


0.0115(3) 


0.0110(3) 


1.0 


10 


0.7019(18) 


14.247(36) 


0.044047(89) 


0.041471(84) 


0.0217(1) 


0.0204(1) 


1.5 


10 


1.0790(12) 


9.268(10) 


0.029886(60) 


0.027989(57) 


0.0348(1) 


0.0326(1) 


3.0 


32 


7.657(35) 


4.179(19) 


0.0014829(25) 


0.0013928(23) 


0.0869(8) 


0.0817(8) 


3.0 


46 


7.114(27) 


6.466(24) 


0.0016151(30) 


0.0015262(28) 


0.0817(6) 


0.0772(6) 


3.2 


60 


9.656(84) 


6.214(54) 


0.0008848(18) 


0.0008400(17) 


0.0825(14) 


0.0783(14) 


3.4 


82 


13.08(17) 


6.267(81) 


0.0004786(16) 


0.0004549(16) 


0.0819(25) 


0.0778(20) 



In Table II we list the run-parameter values used and the main results of simula- 
tions without a topological term. We chose three values of (5 (0.5, 1.0 and 1.5) in the 
strong coupling region, which correspond to correlation lengths of one lattice spacing 
or less. In the weak coupling region, we also chose three (5 values, corresponding to 
correlation lengths of more than seven lattice spacings. Calculations were performed 
on square lattices of size L, chosen such that the condition L/£ > 6 is fulfilled in or- 
der to avoid finite size effects. The condition L/£ > 6 is chosen to be approximately 
the same as that in Ref. 8 ). With this value, a clear scaling plateau of the topological 
susceptibility (as discussed below) has been found. There are finite size effects even 
for this physical size, but the effects are sufficiently small to see the scaling behavior 
for different values of (3 ( (3 = 3.0, 3.2 and 3.4 ), where the effects are expected to be 
similar. The only exception to this is the run performed with (3 = 3.0 and L = 32. 
The difference between the results of this particular run and that on a larger lattice 
clearly demonstrates the existence of finite size effects on smaller lattices, which is in 
accordance with similar observations in Ref. 8 ) . For each run, several million sweeps 
were carried out. 
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During the simulations, we measured the topological charge using two defini- 
tions. One is the standard geometrical charge Qst defined in Eq. (2-2). In addition, 
we measured the FP topological charge Qpp , which is defined as in Ref . 8 ) on a first 
finer level of a multigrid, calculated using a parametrization of the FP field. The 
topological susceptibility is obtained through the relation 

*' - f> , (3,) 

for which numerical values are given in Table II. As in Ref. 8 ), we find that Xfp is 
somewhat smaller than Xst • I n order to study scaling, we consider the dimensionless 
combination x*£ 2 ) whose values are also listed in Table II. We find this quantity to 
be approximately constant for both cases: Xst? 2 ~ 0-082 and Xfp£ 2 ~ 0.077 for the 
three weak couplings. These values are 6 to 12 % larger than the results in Ref. 8 ** but 
display a clear scaling plateau in the same region of £. This confirms that the new 
simpler parametrization of the FP action, even without use of the FP topological 
charge, has the same good scaling properties. This becomes even clearer when one 
compares the behavior with the results for the standard action case, 17 ) where this 
plateau has never been found. Since a combined use of Qpp and the set method 
requires too much computation time, in the simulations discussed below, we used 
Qst for measurements of the topological charge. 

3.2. Scaling behavior of P{Q) and £(Q) 

In this subsection we present the results for P{Q) and £(Q) and discuss the 
scaling behavior of these quantities. 

3.2.1. P(Q) 

Calculations of the topological charge distribution P(Q) using the set method 
were performed for various values of the coupling constant (5 and for various lattice 
volumes V = 1? . An overview is given in Table III, where we also indicate the range 
<2min - Qmax for which P(Q) is calculated in each case. 

Table III. Parameter values used in simulations of CP 3 FP with the set method. 
/3 L: Q 

min ^ max 



0.5 


6 


0-12 


10 


0-18 


















1.5 


8 


0-18 


12 


0-24 


24 


0-39 














3.0 


4 


0-4 


6 


0-12 


8 


0-18 


12 


0- 


-30 


18 


0- 


-54 




24 


0-58 


28 


0-3 


32 


0-94 


34 


0- 


-3 


38 


0- 


-27 




42 


0-3 


46 


0-33 


50 


0-3 


56 


0- 


-3 


96 


0- 


-60 


3.4 


22 


0-18 


44 


0-27 


58 


0-36 















We carried out an extensive study of the lattice size dependence with (3 = 3.0, 
which we have seen in the previous subsection already lies in the scaling region. 
Lattice sizes were chosen in a wide range, beginning from very small ones with 
L/£ < 1 up to large ones with L/£ ~ 13. This study was supplemented by additional 
runs using the larger value of (3 = 3.4. This was done in order to make more 
detailed checks of scaling. Lattice sizes here were chosen so that the values of L/£ 
are approximately the same as those used in the (3 = 3.0 case, resulting in the 
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identifications 22 -£4> 12, 44 24 and 58 44> 32. In order to study the differences 
between the weak coupling region and the strong coupling region, some additional 
simulations were also performed with the strong couplings (3 = 0.5 and 1.5. Since 
for strong coupling the correlation length is too short, the physical lattice size L/£ 
cannot be taken small. For these studies we have therefore been restricted to large 
physical lattice sizes in the range 13 < L/£ < 22 for (3 = 0.5 and 7 < L/£ < 22 for 
/3 = 1.5. 

The statistics of simulations with (3 = 3.0 consisted of the order of 10-20 million 
sweeps per set for the cases in which several sets were employed. Higher statistics 
were obtained, with approximately 50 million sweeps, for simulations with L = 34, 
42, 50 and 56, for which only one set was employed in order to study the behavior 
around Q = 0. The highest statistics, with 150 million sweeps, were obtained for 
L = 46 in the first set, Q = - 3. Increasing the statistics for this case, particularly 
in the first set, was crucial to reduce errors in the free energy density F{9) for 9 
close to 7r. Statistics for (3 = 0.5, 1.5 and 3.4 consisted of about 2 million sweeps 
per set. For error analysis, we employ the jack-knife method by dividing the runs 
into approximately 100 subsets. Since autocorrelation times for Q are at most a 
few of tens of time units, 2 million sweeps were sufficient to make the measurements 
independent and the errors reliably estimated. 




(Q/L) 



Fig. 1. Probability distribution of the topological charge for CP 3 FP with [3 = 3.0 and various 
values of L. Only the results with Q/L < 1 are plotted. The size of error bars is such that they 
cannot be seen at this scale of plotting. The curves connecting the data points are included as 
a guide to the eye. 



In Fig. 1 we plot the measured topological charge distribution P(Q) with (3 = 3.0 
for various lattice sizes L. For convenience and to present data in a compact manner, 
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we normalize the topological charge with the lattice size L. Even a rough inspection 
of the figure reveals that the data do not exhibit Gaussian behavior, which would be 
represented by straight lines with this method of plotting. In particular, the data 
for small lattice sizes exhibit a clear curvature, while they tend to straighten out to 
some extent as the lattice size increases. However, fits with a Gaussian form turn 
out to be extremely poor for all the cases. 

In Fig. 2 we compare the charge distribution obtained with the two weak cou- 
plings (3 = 3.0 and 3.4 for a scaling check of P(Q). For such a comparison to be 
possible, we have now to normalize the charge with the physical lattice size L/£, for 
which we use measured values of the correlation length £. We find that the data at 
the two couplings exhibit the same behavior and lie roughly on the same line. The 
remaining differences are most visible towards larger values of the topological charge. 
This may be explained by the following two reasons. One is that the physical lattice 
sizes can differ slightly. Values of L/£ are 1.687(6) (L=12), 3.37(1) (L=24), 4.50(2) 
(L=32) for (3 = 3.0 and 1.68(2) (L=22), 3.36(4) (L=44), 4.43(6) (L=58) for (3 = 3.4. 
The other is that £ is chosen from Table II for each j3 . We take £ to be 7.114(27) for 
(3 = 3.0 and 13.08(17) for (3 = 3.4. If the measured values of £ were used for each L, 
one would expect at most a 10 % difference among the values of £ in this range of L, 
but the quality of the agreement of the values of P(Q) in the figure does not change. 
In a strict sense, we are not able to obtain rigorous scaling behavior, but the quality 
of the agreement of the two different values of (3 is still remarkably good compared 
to the other two models, CP 1 FP and CP 3 ST. This point is discussed in detail in the 
following section (see Fig. 11). For an argument of more rigorous scaling, we need 
to use measured £ for each L and use Qfp instead of Qst, but here we think the 
results in Fig. 2 are sufficient for the present discussion. 

3.2.2. Correlation length 

Simulations with the set method allowed us to determine the dependence of the 
correlation length on the topological charge. Correlation functions can be measured 
separately for different topological sectors, and correlation lengths £(Q) can be ex- 
tracted from their fall-off at large time separations. The results are shown in Fig. 3, 
where we plot correlation lengths as functions of the charge density. We normal- 
ize the charge density by the correlation length £(9 = 0) obtained from simulations 
without a topological term discussed in Sec. 3.1. 

In the left panel of Fig. 3, we plot the data obtained with (3 = 3.0 and different 
lattice sizes. As expected, the data fall on a universal curve if plotted as functions 
of the topological charge density. We find a decrease of the correlation length with 
increasing topological charge density. The reason for this is that with an increasing 
number of topological objects within the same volume, the configuration has to 
become less ordered. At the point where the topological charge density Q/(L/£(9 = 
0)) 2 is equal to 1, the correlation length £(Q) has dropped to less than half of its value 
at zero topological charge. In the right panel of Fig. 3 we compare representative 
data at two different values of the coupling constant. We see clear scaling behavior of 
£(Q)/L. For smaller values of £(Q), on the other hand, we have observed violations 
of the scaling property. The curve for (3 = 3.0 starts to deviate at Q/(L/£(0 = 
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Fig. 2. Scaling test of the charge probability distribution for CP 3 FP at two values of the coupling 
and several volumes. 



0)) 2 ~ 0.7 where £ ~ 4. This, in turn, indicates that even a use of the FP action is 
unable to suppress the lattice defects for £ < 4, which is in agreement with similar 
previous observations. 8 )' 20 ) 

From these scaling checks, we conclude that the observed behaviors of P(Q) 
and £(Q) for (3 = 3.0 represent continuum properties and are not caused by lattice 
artifacts. These are compared with the results for CP 1 FP and CP 3 ST actions in the 
following section. 

3.3. Free energy density and expectation value of the topological charge 

In Fig. 4 we show results for the free energy density obtained with Eq. (2T3) for 
(3 = 3.0 and several lattice sizes L. Central values are indicated by solid lines, while 
the one-sigma error bands, determined with a jack-knife analysis, are represented by 
dashed lines. We obtain smooth curves, and no "flattening" is observed, as in some 
other works. n )' 18 )' 19 ) The only exception is the simulation with L = 96 (not shown 
in the figure), where the calculation of F{9) breaks down at 9 w 1, because the 
partition function turns out to be negative (while still consistent with zero within 
the error bars). Before this breakdown, the error bars become very large. 

Note that the physical free energy densities for the corresponding lattice sizes of 
[3 = 3.0 and (3 = 3.4 fall on the same curve, which is a consequence of the scaling of 
P(Q) shown in Fig. 2. We observe a strong dependence of the free energy density 
on the lattice size. F{6) increases with increasing L and for 9 < 1, already, nearly 
asymptotic values are obtained for our largest lattice sizes. This is, however, not yet 
the case for 9 close to it. This has the consequence that curves become steeper as 
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Fig. 3. Fixed charge sector correlation length as a function of the physical charge density for CP FP. 
In the left panel we compare various lattice sizes at j3 = 3.0 and in the right panel we show 
scaling between (3 = 3.0 and 3.4. 
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Fig. 4. Free energy density F{8) for CP 3 FP for [3 = 3.0 and several lattice sizes L. 
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the volume increases. Does this mean that F{9) develops a peak at 9 = w, which 
would be an indication of a phase transition? 

To answer this question, we consider the expectation value of the topological 
charge obtained with Eq. (2-14). We show results obtained for (3 = 3.0 in Fig. 5. 
As seen there, {Q)e vanishes both for 9 = and n. In between these values, it has 
a peak, which we observe to move slowly away from 9 = ir/2, obtained at small 
lattices, towards it. Fig ure 6 displays the volume dependence of 9 mStX > defined to 
be the position of the peak of (Q)g. This clearly shows how # max moves away from 
7r/2. It is, however, still far from ir even with our largest lattice, and it cannot be 
conclusively determined from our data where the infinite volume limit of 9 max is. 
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Fig. 5. Expectation value of the topological charge (Q)e for CP 3 FP with (3 = 3.0 and several lattice 
sizes L. Results increase from L — 4 (data points at the bottom) to L — 46 (data points at the 
top). 



3.4. Effective power j e ff of P(Q) 

We are now in a position to investigate the behavior of the charge distribution 
in more detail without having to be concerned about discretization effects. This is 
important for obtaining information about the underlying physics. Having observed 
that P(Q) cannot be fit with a simple Gaussian, we turn to a closer examination 
of its local properties. To this end, we calculate the effective power = ^effiQ) 
defined by assuming that P(Q) behaves at the three adjacent charges Q, Q + 1 and 
Q + 2 as the function 

flocai(Q) = Aexp(-BQ^). (3-2) 

Since three values of the charge distribution are used as inputs and the function 
in Eq. (3-2) has three free parameters, the latter can be calculated algebraically 
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Fig. 6. Angle # max of the maximal topological charge for CP FP as a function of the physical lattice 
volume. 



without involving a fit. Trivially, the Gaussian distribution takes the value jejj = 2, 
independent of Q. 

Results for j e ff{Q) with (3 = 3.0 are displayed in Fig. 7 as a function of the charge 
density or charge filling fraction Q/V of the lattice. We plot the same data twice, 
first showing all results up to large filling fractions in the left panel, and then zooming 
in to the region of small filling fractions in the right panel, where data for small (4, 6 
and 8) and large (96) lattice sizes are omitted for clearer visibility. We observe that 
for almost all the combinations of lattice size L and topological charge Q, the values 
of ~f e ff differ considerably from 2. This explains why we could not fit P(Q) with a 
Gaussian. In fact, P(Q) is not even well described locally by a Gaussian for most 
of the cases. We find, however, that the values of 'YeffiQ) are bounded between the 
two limits 1 and 2. The values found at Q = depend on the lattice size. When the 
topological charge is increased, we find that results from different lattice sizes seem 
to approach a common value of ~ 1.3. This common value is reached, at the 
latest, at a filling fraction of around 1/10. An interesting point to observe is that the 
approach to this common value is from two different sides, depending on the lattice 
size. Data for L = 12, for example, start with jejj = 1.13 at Q = and increase 
with increasing charge. As an example of different behavior, data for L = 38 start 
with 7 ej ^ = 1.88 at Q = and decrease with increasing charge. The lattice size at 
which increasing behavior changes to decreasing behavior is approximately L = 24 
(which corresponds to L/£ ~ 3.5). 

If we increase the topological charge to filling fractions larger than about 1/10, 
the values of j e ff move away from 1.3 and increase in a universal manner. This can 
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Fig. 7. Effective power of the probability distribution for CP 3 FP with f3 — 3.0 and various lattice 
sizes. The left panel shows all results up to a charge filling fraction of 1/4, while the right panel 
is an enlargement of the region with small filling fraction. 



be easily understood to be an effect of approaching the limit of the maximal number 
of topological objects (instantons) that can be placed on a lattice. The probability 
of such densely packed lattice configurations is decreased, which, in turn, leads to 
an enhancement of "Jeff- This point will be investigated in section 5 by studying 
analytical models. 

In Fig. 8 we compare the results for 7 c fj at the two weak couplings, (3 = 3.0 and 
3.4. In order to make this comparison, we have to plot the data as a function of 
the physical charge density instead of the filling fraction. By doing this, the data 
at the two couplings collapse to a common line, showing scaling also for Jeff- This 
finding also confirms the approach from different sides to the asymptotic value for 
different physical lattice sizes. Let us finally focus on the value of j e ff obtained 
at the origin, Q = 0. As observed in Fig. 7, it takes values between 1 and 2 and 
exhibits strong dependence on the lattice size. This lattice size dependence is shown 
in Fig. 9. We clearly observe a crossover from values close to 1 for small lattices to 
values approaching 2 in the limit of large lattices. In this way, we recover a Gaussian 
charge distribution at weak coupling. It is realized, however, only in the limit of 
large volume and only near vanishing topological charge. The effect observed at 
weak coupling is also observed for strong coupling, as shown in Fig. 10. Here we 
find "f e jj to be around 2, in accordance with the observation of Gaussian behavior 
at strong coupling in Refs. 10 ) and 11 ). However, if the filling fraction becomes large, 
we again see an increase of 7 e jg F , so that even values larger than 2 can be obtained. 

It is interesting to note in connection to this behavior that numerical simulations 
in the strong coupling region, where a Gaussian distribution was observed in previous 
studies, 10 )' n ) are confined at the same time in a region of large volumes and in a 
region of small physical charge density. The reason for this is that for a strong 
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Fig. 9. Effective power y e ff at the origin as a function of the lattice size L for CP 3 FP at j3 = 3.0. 

coupling, the correlation length £ is small. The lattice size L, however, should not 
be chosen to be small, so that the lattice can contain more than just a few instantons. 
As a result, the ratio L/£ becomes large, which means that we are in a region of 
large physical volumes. Moreover, even though in the above-mentioned previous 
calculations the range of Q was taken up to fairly large values, they still correspond 
to a small physical charge density £ 2 Q/V , because values of L/£ are large. 

§4. Numerical results for CP^P and CP 3 ST 



It is worthwhile to investigate whether scaling behavior is observed, and if not, to 
what extent it breaks down for other cases of the CP^" 1 model. For comparison, we 
considered a model with CP 1 FP action. As a parametrization of the CP 1 FP action, 
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Fig. 10. Jeff for CP 3 FP at strong couplings. 

we used the parametrization listed in Table 4 of the original paper by Hasenfratz 
and Niedermayer, 5 ) in which 24 local coupling constants are employed. For updating 
configurations, we used the Metropolis algorithm. Typically, we performed one to 
several million sweeps per set, depending on the coupling constant (3 and L. For the 
set method described in section 2, we always chose 4 bins for each set and the trial 
function to be Gaussian. The values of (3 ranged from 0.7 to 1.1 and the lattice sizes 
from 12 to 62. 

In order to study the scaling behavior, we compare the results for pairs of pa- 
rameter doublets (/?, L), as done in the previous section, so that each one of the pair 
is chosen to have approximately the same physical volume, (L/£(# = 0)) 2 . Figures 
11 and 12 show P(Q) and F(9), respectively, for (/?, L)=(0.87,22) and (1.02, 44), 
where L/£{9 = 0) is measured to be 3.2. These data exhibit a large difference, in 
contrast to the CP 3 FP case depicted in Fig. 2, and the scaling is found to be clearly 
broken. In Fig. 13 we compare the correlation lengths £/L as a function of the physi- 
cal charge density Q/(L/£(9 = 0)) 2 . The two sets of data for £/L also exhibit strong 
scaling violation, and large deviation(« 30%) is found for Q/(L/£(9 = 0)) 2 « 0.7. It 
is noted that this behavior is also seen for other choices of pairs with different values 
of L/£(9 = 0), e.g., ((3,L) =(0.7, 12) and (0.9,30), with L/£(d = 0) being 3.6. 

We have found for the CP 1 model that scaling is strongly violated, although 
the FP action is used. This is expected, since it is well known that dislocations do 
harm seriously in the CP 1 model. In our calculations, we have used only the FP 
action, and the additional use of a FP charge might be expected to produce better 
results. However, we suspect that even if we had used the FP charge, the scaling 
nature would not have been improved. This comes from the observation that even 
after the dislocations are eliminated by a combined use of the FP charge and the 
FP action, the scaling behavior of the lattice topological susceptibility xt is strongly 
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violated. 6 )' 7 ) 

We have also calculated the CP 3 standard action for different choices of f3 and 
L. We varied lattice sizes from 18 to 54 and (5 from 1.8 to 2.3 and looked for pairs of 
doublets (/?, L) such that L/£ for each pair is approximately the same. Values of L/£ 
were chosen to be from about 2.7 up to 4.2, which are similar to the values used for 
CP 1 FP and CP 3 FP. For updating configurations, we used the Metropolis algorithm, 
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and the number of sweeps used for each set was approximately same as for CP FP. 
Use of the standard action for the CP 3 model exhibits strong violation, similar to 
the case CP 1 FP, for P(Q) and F(9), as well as the correlation length £/L. This is 
easily understood, since the standard action is very poor in representing continuum 
physics. The behavior of £/L exhibits somewhat small deviation compared to that 
for CP^P. 

§5. Debye-Hiickel approximation 

In Refs. 10 ) and n ) it has been shown that a first-order phase transition exists 
at 6 = it when the topological charge distribution P(Q) is Gaussian, and its volume 
dependence is given by P(Q) ~ exp[— aQ 2 /V], where a is a constant that depends 
on the coupling constant (3. Such a Gaussian charge distribution has been found 
in the region of very strong coupling. As [3 becomes larger for some fixed volume, 
P(Q) has been observed to deviate from the Gaussian form, and thus the first-order 
phase transition gradually disappears. We note that P(Q) for a system that consists 
of instantons and anti-instantons obeying the Poisson distribution also behaves like 
a Gaussian for values of the parameter corresponding to the strong coupling region. 
Therefore, we expect instantons and anti-instantons in such a system to behave like 
a dilute gas in the strong coupling region. It is thus of interest to determine the 
nature of the dynamics of instantons displayed by a system for which P{Q) is not 
Gaussian. In order to investigate this issue, we use the Debye-Hiickel model, 12 ) 
which is based on an instanton quark picture, and in which correlations between 
particles (instanton quarks) are weak. 
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5.1. Instanton quark picture and Debye-Hiickel model 

In this subsection we explain the concept of the instanton quark picture and the 
Debye-Hiickel model (D-H model). In order for this paper to be self-contained, we 
also give a summarized overview of the results obtained in previous works. 12 )' 13 ) In 
Ref. 13 ) , Fateev, Frolov and Schwarz analyzed Euclidean Green's functions and the 
partition function to investigate how instantons affect the dynamics of the CP W_1 
model. The partition function of this model in the continuum is defined by 



(5-1) 



Z = f Vzml[S(J2 \*a{x)\ 2 - l)e~ s , 
and the action is defined by 

N 

S = (3 d 2 x ]T{| ^z a (x)| 2 + (z a {x)d^z a (x)) 2 }, (5-2) 

J a=l 

where ~z a denotes the complex conjugate of z a . The (N — l)-dimensional complex 
projective space is defined by introducing a field as 

MO = MO/MO; a = 1, 2, ■ • • , N - 1 , where z N (() + 0, (5-3) 
in the complex plane ( = x + iy. As a next step, the action is rewritten in the form 

S = (3 J d x — 2 2^ \dnM 



where p 2 = J2 a =i \ w a\ 2 and d 
by 



(54) 

d(d(. A general g-instanton solution is given 

q 



«(C)=CaII(C-«a), 



(5-5) 



where c a and a l a are complex parameters. The superscript i (i = 1,2,- ■• ,q) labels 
instantons, and the subscript a ( a = 1, 2, • • • , N) indicates the degrees of freedom of 
the field z a . The authors of Ref. 13 ) investigated how the system behaves when the 
field fluctuates around the g-instanton solution 21 )' 22 ) and showed that the partition 
function takes the form 

. N q N 

Z = Y.W q Wr N UU d2c - d2a a^{-n q (a,c)/T}6{j2\^\ 2 -l), (5-6) 

q J a =li=l a=l 

where W is a constant that depends on the topological charge, and T = 1. H q in 
Eq. (5-6) is given by 



H q {a,c) = ^- f d 2 x\\ogpd 2 \ogp-\og(j{\c a \ 2 «X{\a i a -a^\ 

^ L a i>j 



(5-7) 



Equation (5-6) can be interpreted as the partition function of a system of two- 
dimensional classical particles with interaction energy given by Eq. (5-7) in the grand 
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canonical ensemble. The classical particles are at positions a l a and interact with each 
other. Note that T plays the role of a temperature. 

In the particular case N = 2, i.e. for the 0(3) non-linear sigma model, the grand 
partition function is given by 

Z' = J2W«(q\)- 2 f Hd 2 a\d 2 alexp{-n' q (a)/T}, (5-8) 
q J i 

H' q {a) = - ]T log \a\ - a{ | 2 - ]T log \a\ - a\f + ]T log \a\ - a\f . (5-9) 

i<j i<j i,j 

We see that in the Hamiltonian (5-9), a particle possesses a "charge" (a = 1,2) 
and that particles of equal charge interact repulsively, while those of opposite charge 
interact attractively. Thus we can interpret this model as a two-dimensional classical 
Coulomb system that consists of 2q particles with positive and negative charges. 
Furthermore, if the locations {a l a ;i = 1,2, ■■■,q, a = 1,2} satisfy the conditions 
\a\ — a\\ <C \a\ — a{\, \a\ — a J 2 \(i ^ j), we can regard \{a\ + a % 2 ) as the position and 
\\a\ — 0-2 1 as the size of the i-th instanton. Due to the interaction (5-9), a pair of 
particles with opposite charges tends to make up an instanton with neutral charge. 
For this reason, the particles are called "instanton quarks" . 13 ) 

In the case of a g-anti-instanton configuration, the solution is given by 

Za(C) = c a f[(C-V a ), a = 1,2, 
i=i 

and the partition function has the same form as in Eq. (5-8) but with the variables 
a l a replaced by b l a . In analogy with the 0(3) non-linear sigma model, the CP^ -1 
model can be interpreted as a classical system that consists of instanton quarks 
with "multicolors" (a = 1, 2, ■ ■ ■ , N) such that instanton quarks with the same color 
interact repulsively and those with different colors interact attractively. 

Recently Diakonov and Maul 12 ) investigated the effects of instantons and anti- 
instantons for the CP^ -1 model by using the instanton quark picture. They proposed 
a configuration with N + instantons and 7V_ anti-instantons by assuming a product 
ansatz of the form 23 ) 

N+ N- 

*a(o=n(c-««)iiG-£). ( 5 - iq ) 

1=1 j=l 

where a = 1,2,---,N. Since Eq. (5-7) cannot be analytically solved for N > 2, 
they simplified the multi-body interaction into an interaction of the two-body type. 
Furthermore, they applied the Debye-Huckel approximation to this model; i.e., they 
assumed that the correlations between instanton quarks are very weak and individual 
quarks interact with a mean field. The partition function turns out to be 



Z{9)= ]T Z(N + ,N_)exp{i6(N + -N_)}, (5-11) 
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where the partition function Z(N + , N-) is given by 



Z(N + ,N_) =exp 



N 2 V d 
~4(N - 1) 



0< 



n + + n_ 



log 



A 2 



--y/(n + -n. 



) 2 + 4n+n_/3' 2 log — 



(V d A 2 ) N{ - N + +N -) 
{N + \) N (N_\) N 



+ 1 



(5-12) 



Here 



n± 

K± 

A 



1/T 

N 



N ± 



N-lV d 

27r/3[(n+ + n_) ± ^/(n + - n_) 2 + 4n+n^/3' 2 
27le 7B . 



In the above, /?' is a coupling constant between instanton quarks and anti-instanton 
quarks introduced by the ansatz represented by Eq. (5-10), and Vd(= L d ) is the 
volume of the system in the continuum. A is a parameter like Aqqd, an d is the 
Euler number. Note that T plays the role of a temperature ( here T = 1). 
We calculate P(Q) from Eq. (5-12) as 

Z{9)= ]T Z{N + ,NJ)e X p{W{N + -NJ)} 

= E E S N+ - N _, Q Z(N + ,N-)exp{ie(N + -N-)} 

Q N + ,N_ 

= ^exp{^Q}P(g), (5-13) 
Q 

where 

P(Q)= E Sn + -n^,qZ(N + ,N_). (5-14) 

We carefully consider the behavior of P(Q) in the next subsection. 

5.2. Analysis of P(Q) in terms of the instanton picture based on instanton quarks 

We now discuss the topological charge distribution P(Q) in terms of instantons 
(anti-instantons) . Before referring to P(Q), let us first consider the behavior of 
Z(N+, N-). In Fig. 14 we give a two-dimensional plot of Z(N+,N-) obtained from 
Eq. (5-12) for L d = 4, (3 = 1.0 and f3' = 0.0. These values of the parameters were 
chosen to be the same as those in Ref. 12 \ From this point, we set A to unity. 12 ) 
Figure 14 displays typical behavior of Z(N + , iV_). This function is symmetric under 
interchange of N+ and iV_ and has a maximum at some non-zero value of iV_|_(= 
N-). The region around the maximum gives a dominant contribution to the charge 
distribution P(Q) at Q = 0. It is interesting to see that a fairly large number of 



22 



R. Burkhalter, M. Imachi, Y. Shinno and H. Yoneyama 




~i 1 1 r 

10 20 30 40 50 60 70 80 

Nm 



Fig. 14. Two-dimensional plot of the partition function log Z(N + , N- ) for N=4, Ld — 4, j3 = 1.0 
and = 0.0. Only results with 1 < N p ,N m < 80 are plotted. Here N p and N m are the number 
of instantons and anti-instantons, respectively. 



instantons (N + = N- ~ 30) contribute to P(0). In Fig. 15 we compare Z(N + , N-) 
for two different values of (3' {(3' = 0.5 and (3' = 0.0) in the case Q = (N+ = iV_). 
We find that Z(N + , iV_) has only weak (3' dependence. In fact, this is true for other 
choices of the parameters as well. This result shows that the correlations between 
instanton quarks and anti-instanton quarks can be ignored. 12 ) For this reason, we 
set = 0. Figure 16 shows how Z(N+, N-) depends on /3 for a fixed volume. We see 
that as (3 becomes larger, the position of the maximum of Z(N+,N-) shifts to the 
right and configurations with larger numbers of instantons and anti-instantons are 
generated. Note that (3 = 1.0 corresponds to a physical system. For fixed (3, on the 
other hand, the peak of Z{N + ,N^) shifts upward and to the right with increasing 
volume. This is intuitively understandable. 

We are now in a position to study P{Q). In Fig. 17 we plot P(Q) obtained from 
Z(N + , iV_) for (3 = 1.0 and = 4. At first sight, P{Q) looks like a Gaussian. In 
order to investigate the behavior of P(Q) in more detail, we calculate the effective 
power -f e ff introduced in subsection 3.4. We recall that in Monte-Carlo simulations, 
j e ff is nearly equal to 2.0 at Q/V = and decreases slowly with increasing charge 
for large volumes, while for small volumes, j e jj is near unity at Q/V = and 
increases rapidly with increasing charge. It is interesting to investigate the cause of 
this difference and whether it can be accounted for by the D-H model. In Fig. 18 
we display 7^ as a function of the physical charge density Q/Vd, where 7^ is 

obtained in the D-H approximation. Obviously 7^' is different from 2.0 (Gaussian). 
Even if we take (3 1 7^ 0, the results do not change significantly. Compared to Fig. 7, 
the D-H approximation is found to reproduce the behavior of j e g for large volumes. 

Let us now turn to small volumes, for which we believe that finite size effects 
are relevant for P(Q) on a lattice. Since instantons with a size smaller than the 
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Fig. 15. The logarithm of Z(N+, JV_) for the case Q=0 with j3 = 1.0 and L d = 4. The results for 
= 0.0 are represented by solid dots and the results for /3' = 0.5 by crosses. 
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Fig. 16. The logarithm of Z(N+,N-) in the case Q=0 with Ld = 4 and /?' = 0.0. (3 is chosen as 
0.1(«), 0.5(+) and 1.0(x). 



lattice spacing fall through the lattice, the number of instantons is bounded from 
above (= Nm'ax) in finite volumes. Contrastingly, the D-H model is defined in the 
continuum, and for it, an infinite number of instantons can be placed in finite volume. 
In order to reproduce the lattice results with this model, it is necessary to restrict 
the upper limit (= N max ) to a finite value in the summation J2n+,N- N-) in 
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Fig. 17. The logarithm ol P(Q) as a function of Q 2 in the Debye-Hiickel approximation for (3 = 1.0, 
/?' = 0.0 and L d = 4. 
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Fig. 18. 7^ obtained from P(Q) for (3 = 1.0, /?' = 0.0 and L d = 4. 



Eq. (5-14). Let us assume that N max niiniics N max . In Fig. 19 we display the results 
for = 1.0 and = 4 with three different choices of N max . For large iV maa; (> 40), 
starting at Q = 0, 7^ gradually decreases from 2.0 and starts to diverge at some 

non-zero value of Q/Vd- For smaller values of N rnax (= 25), 7^ is close to unity at 
Q = and increases rapidly as a function of Q. At some small value of the filling 
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25 



gamma 




Q/v d 



Fig. 19. Same as in Fig. 18, but with an upper limit N„ 
N max is taken to be 25(»), 50(x) and 75(+). 



in the summation ^2 N+ N Z(N+ 



fraction, y^' overshoots 2.0. The behavior seen in Fig. 19 is quite similar to that 

seen in Fig. 7. It is noted that 7^ at Q/Vd = gradually shifts up from 1.0 up 
to 2.0 if we increase N max for small values, in accordance with the behavior seen in 
Fig. 9. 

We find by varying N max and by taking various values of (3 and L that when 
Nmax is chosen to be approximately equal to or less than the location of the peak 
of the distribution Z(N + ,N_) (see Fig. 14), 7^ at Q = starts at a value nearly 

equal to 1, and when N rnax is larger than the location of the peak, 7^ starts at 

2. For the latter case, 7^ slowly decreases as a function of Q/V^, and diverges 

at some value of Q/Vd- As N max becomes even larger, 7^ begins to diverge at a 
larger density Q/Vd- These observations allow us to conclude that the behavior of 
P(Q) found in Monte Carlo simulations reflects the dynamics of weakly correlated 
instanton quarks. 

5.3. Analysis of P(Q) in terms of the Poisson distribution 

Before ending this section, we study j e jj obtained from a system that consists of 
instantons and anti-instantons obeying the Poisson distribution. This is equivalent 
to a dilute gas system (DGA) of instantons and anti-instantons for which P(Q) is 
given by 

P (Q)= E 2 N+ + N. (N ^~1 )I PN + +N - 6N +- N -' Q > (5 ' 15) 
N+.N- + ' 
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P = —e~ x 
n! 

where A/2 is the average number of instantons (anti-instantons) . 
Equation (5-15) leads in the limit A — > oo to 

P(Q)~exp{-^Q 2 }. (5-16) 

With the identification A = X, Eq. (5-16) can be regarded as the distribution in the 
strong coupling region of MC simulations, where a is a constant that depends on (3 
appearing in P(Q) ~ exp[— yQ 2 ] in such a manner that a — > as — > 0. 

(p) 

Let us study how ^ e g behaves for A < oo. Figure 20 displays "f e ff(= 7^ ) as a 
gamma 
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+ 
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•:N max =20 




+:N max =40 




x:N max =80 



Q/lambda 



Fig. 20. j^jp f° r ^ = 40, where A is considered to be proportional to Vd- 

function of Q/ A for A = 40. We find that the behavior of 7^ is very similar to that 

of 7g«^ for various values of N max , where N rnax is the upper limit of the summation in 
Eq. (5- 15). Recalling that in the D-H model, correlations between instanton quarks 
are weak, while in the DGA there are no correlations between instantons, we can 
conclude that the D-H model also has little correlation between instantons. 

§6. Conclusions and discussion 

(1) We have studied CP^ -1 models with a topological term. In order to obtain a 
description of continuum physics, we have employed an FP action and investi- 
gated scaling properties of quantities such as P(Q), F(6) and the correlation 
length £(Q)/L. For CP 3 FP we have observed good scaling behavior, while for 
CP 1 FP and CP 3 ST we have found strong violations of scaling. 
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We have investigated the effective power ^ e jj of P(Q) for CP 3 FP. At a fixed 
value of (3, j e ff increases from 1.0 for small lattices as Q increases, while j e ff 
decreases from 2.0 for large L. The values of for both large and small 
lattices reach a common line at some small value of the filling fraction. When 
finite size effects become significant, the value of ^y e jj starts increasing beyond 
2.0. 

We have studied the behavior of for the analytical models, the Debye- 
Hiickel approximation to an instanton quark gas of the CP^ -1 model, and the 
Poisson distribution of an instanton gas. We have found that in these cases, 
exhibits the same behavior as that obtained in Monte Carlo simulations. Finite 
size effects, which emerge by packing instantons into a finite volume, show up 
as an increase beyond 2.0 of "Jeff- These observations allow us to conclude 
that P(Q) obtained in MC simulations describes the dynamics of very weakly 
correlated instantons. 

Gaussian behavior of P(Q) seems to be realized when volumes are large and the 
physical charge density is small. In the strong coupling region, these conditions 
are easily satisfied, because correlation lengths are very short in this region. 
As a consequence, there exists a first-order phase transition at 9 = it. In the 
weak coupling region, on the other hand, "f e ff tends to be 2.0 only at vanishing 
topological charge when the volume increases. The expectation value (Q)q 
develops a peak, and there is the tendency that the peak moves away from tt/2 
towards it as the volume increases. However, we could not obtain conclusive 
results about existence of a phase transition in the infinite volume limit, because 
the above conditions are difficult to satisfy. 

When the volume is small, the value of j e jj increases from 1.0 as a function 
of the charge in all the cases we have investigated, not only for simulations 
but also for analytical models. It is an interesting question why is always 
bounded from below by 1.0. We think that perhaps this would be associated 
with some fundamental property of probability theory. 

Although the analytical models can explain the behavior of qualitatively 
well, 'Jeff for small volumes increases more rapidly than in MC simulations. This 
behavior might be due to the sharp cut-off (N max ) used in the summation in 
Eqs. (5-14) and (5-15). We have tentatively smeared the cut-off by introducing 
a Gaussian function. As a result, we have observed that the increase of j e jj 
becomes somewhat less rapid, and it tends to follow the common line before 
diverging. However, in order to draw a definite conclusion, a more systematic 
study is needed. 

The fate of the first-order phase transition at 6 = tt in the strong coupling 
region is a relevant issue. As discussed in subsection 3.3, calculations of the 
free energy and the expectation value (Q)g have huge errors already at L = 96, 
and thus we cannot address this issue. To increase lattice size while keeping 
errors reasonably small requires an exponentially increasing number of sweeps. 
We may need some novel algorithm to overcome this problem. 
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